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Abstract 

^v^j . Complementary measurements have been performed for the level den- 

■ sity and 7 strength function in 162 Dy using inelastic 3 He scattering. Com- 

f^*) [ paring these results to previous measurements using the 163 Dy( 3 He,a) re- 

^SJ ■ action, reveals that the measured quantities above 1.5 MeV do not depend 

significantly on the nuclear reaction chosen. 
PACS number(s): 21. 10. Ma, 21.10.Pc, 25.55.-e, 27.70. +q 

G\ 

1) ; 1 Introduction 

^) ' Nuclear level densities have recently gained new interest. When earlier studies 

of level densities were mainly based on counting levels close to the ground state 
and neutron resonance spacing at the neutron binding energy [Q, , a variety 
of new methods and experimental results are available today. A more recent 
compilation of all existing data on level densities §1 includes level spacing data 
of several other reactions involving light particles up to A — 4 as well as results 
from Ericson fluctuation measurements. Recently, experimental level densities 
in 69 As and 70 Ge over a large excitation energy interval of 5-24 MeV have been 
reported H, obtained from proton evaporation spectra of 12 C induced reactions. 
Also the Oslo cyclotron group has reported on a new method to extract level 
density and 7 strength function from primary 7 spectra (see Q for the basic 
assumptions and || for the method). This method has the advantage that 
the level density is deduced from 7 transitions, thus the nucleus is likely to be 
thermalized and the measured level density is supposed to be independent of 
the formation mechanism of the excited nucleus. Several applications of the 
method are reported in |7|, ||, ||, [l^] . 

The experimental progress has been accompanied by new theoretical devel- 
opments, with respect to the first analytical nuclear level density formula pro- 
posed by Bcthe IllT . Level densities have been studied for finite temperatures 



1 



within the BCS model 12, [Tj|. Today, Monte Carlo shell model calculations 
£3 Ell are a ^ e ^° estimate nuclear level densities |l6| for heavy mid shell nuclei 
like ""Dy jl?]]. Also more schematic approaches like binomial level densities [ [l8[ 
have been revived lately. Important applications of the theoretical and experi- 
mental efforts are calculations of the nucleon synthesis in stars, where the level 
densities are inputs in large computer codes and thousands of cross sections are 
estimated |l9| . 

Also the present knowledge of the 7 strength function is poor. Although the 
strengths can be roughly calculated by the Weisskopf estimate, which is based 
on single particle transitions (see e.g. |2(]]), some transitions deviate many orders 
of magnitude from this approximation. A compilation of average 7 transition 
strengths for dipole and electric quadrupole transitions can be found in pj| . 
The uncertainty of the 7 strength function concerns the absolute value and 
the 7 energy dependence. For El transitions one assumes that the 7 energy 
dependence follows the Giant Dipole Resonance (GDR) (7,7') cross section. 
This is, however, to be proven. 

In this work, we determine the level density and the 7 strength function for 
162 Dy for energies close up to the neutron binding energy B n . By comparing the 
present data, which were obtained from the 162 Dy( 3 He, 3 He'7) 162 Dy reaction, 
to previous data |22|, ||, which were obtained from the 163 Dy( 3 He,a7) 162 Dy 
reaction, we can test if the basic assumption of our analysis method is fulfilled. 

This main assumption is that the 7 decay pattern from any excitation energy 
bin is independent of the population mechanism of states within this bin, e.g. 
direct population by a nuclear reaction, or indirect population by a nuclear 
reaction followed by one or several 7 rays. Since the 7 decay probabilities of an 
excited state are independent of the populating reaction, the assumption above 
is generally equivalent to the assumption that the same states are populated 
equally by the direct and indirect population mechanisms. One can now imagine 
several cases where this assumption might be invalid. 

Firstly, thermalization time might compete with the half life of excited states, 
and the selectivity of the direct population by a nuclear reaction will be reflected 
by a different 7 decay pattern with few and relatively strong 7 transitions com- 
pared to a statistical spectrum which is the expected 7 decay pattern after 
complete thermalization. 

Secondly, direct population might populate states with different exact or 
approximate quantum numbers like spin or parity than indirect population. 
Since states with different exact or approximate quantum numbers do not mix 
at all or very weakly in the latter case, the ensemble of populated states after 
thermalization will differ for the two population mechanisms and therefore one 
can expect different 7 decay patterns. 

It is very difficult to judge where the assumption of the method is appli- 
cable and how good this approximation is. Below, we will, by comparing two 
different direct population mechanisms represented by two different nuclear re- 
actions, investigate in which excitation energy interval the assumption might 
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break down. 



2 Experiment and data analysis 

The experiment was carried out at the Oslo Cyclotron Laboratory (OCL) us- 
ing the MC35 Scanditronix cyclotron. The beam current was ~1 nA of 3 He 
particles with an energy of 45 MeV. The experiment was running for a total 
of 2 weeks. The target was an isotopically enriched 95% 162 Dy self support- 
ing metal foil with a thickness of 1.4 mg/cm 2 glued on an aluminum frame. 
Particle identification and energy measurements were performed by a ring of 
8 Si(Li) telescopes at 45° relative to the beam axis. The telescopes consist of 
a front and end detector with thicknesses of some 150 /zm and 3000 /im re- 
spectively, which is enough to effectively stop the ejectiles of the reaction. The 
7 rays were detected by a ball of 27 5"x5" Nal(Tl) detectors (CACTUS) @ 
covering a solid angle of ~15% of 47T. Three 60% Ge(HP) detectors were used 
to monitor the selectivity of the reaction and the entrance spin distribution of 
the product nucleus. During the experiment we collected besides data for the 
162 Dy( 3 He, 3 He') 162 Dy reaction, where results are presented in this work, also 
data for the 162 Dy( 3 He,a) 161 Dy reaction, where some results were presented in 
[|| [l0|. A comprehensive description of the 163 Dy( 3 He,a7) 162 Dy experiment, 
which we will compare our findings to, can be found in p4| . 

In the first step of the data analysis, the measured ejectile energy is trans- 
formed into excitation energy of the product nucleus. In Fig. [l] the raw data 
are shown. In the next step, the 7 spectra are unfolded for every excitation 
energy bin using measured response functions of the CACTUS detector array 
[ p5| . In Fig. H the unfolded data are shown. In the third step, the primary 7 
spectra for every excitation energy bin are extracted from the unfolded data 
by the subtraction technique of Ref . |p6| . In Fig. || the primary 7 spectra are 
shown. 

In the fourth step, we extract level density and 7 strength function from the 
primary 7 spectra. The main assumption behind this method is the Axel Brink 
hypothesis |27| p8| 

r(E x ,Ej<xF(Ej e (E f ) (1) 

with Ef = E x — Ey. It says that the 7 decay probability in the continuum 
energy region represented by the primary 7 spectrum T is proportional to the 
level density g and a 7 energy dependent factor F. The level density and the 7 
energy dependent factor are estimated by a least % 2 fit to the experimental data 
[pj . In Fig. U the experimental data including estimated errors |J are compared 
to the fit according to Eq. (|l|). 

The data are fitted very well by the theoretical expression of Eq. (Q) . This 
is a remarkable example for the validity of the Axel Brink hypothesis. However, 
it can never completely be ruled out, that a minor portion of the primary 7 
matrix cannot be factorized into a level density and a 7 energy dependent factor. 
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One might also encounter large fluctuations in these quantities at very low 
level densities around the ground state or when considering highly collective 7 
transitions and single particle 7 transitions at similar 7 energies. 

Since the least \ 2 fit according to Eq. ([!]) yields an infinitely large number 
of equally good solutions, which can be obtained by transforming one arbitrary 
solution by 

g(E x ~E 7 ) = g(E x - E 7 ) Aexp(a[E x - E^]) (2) 
F(E y ) = F(E^)Bexp(aE y ), (3) 

[|| we have to determine the three parameters A, B and a of the transformation 
by comparing the results to other experimental data. We fix the parameters A 
and a by comparing the extracted level density curve to the number of known 
levels per excitation energy bin around the ground state |^9| and to the level 
density at the neutron binding energy B n calculated from neutron resonance 
spacing data |3(|. Since the procedure is described in detail in Ref. ||, we 
only show in Fig. |5| how the extracted level density curve compares to other 
experimental data. 

The parameter B could now in principle be fixed by comparing the extracted 
7 energy dependent factor F to other experimental data of the 7 strength func- 
tion. However since data are very sparse and the absolute normalization of 7 
strength function data is very uncertain, we give the 7 energy dependent factor 
in arbitrary units. 

3 Results and discussion 
3.1 The level density 

We compare extracted level densities of 162 Dy from two reactions, namely 
162 Dy( 3 He, 3 He'7) 162 Dy and 163 Dy( 3 He,a7) 162 Dy. While level densities from 
the latter reaction were already published in ^2], Q using approximate ex- 
traction methods, and in Q in the present form, data from the first reaction 
are shown here for the first time. Figure |^ shows the relative level densities, 
which are calculated by dividing the extracted level densities by an exponential 
Cexp(£;/T) with T = 580 keV and C = 10 MeV -1 in our case. One can see that 
both level densities agree very well within 10% in the excitation energy interval 
1.5 MeV to 6.5 MeV. This result is very encouraging, since level densities are 
generally only known within an error of ±50-100%. Above 6.5 MeV the errors 
are too large in order to make conclusive observations. Below ^1.5 MeV the two 
level densities differ dramatically from each other. In Fig. || one can see that the 
extracted level density from the 163 Dy( 3 He,a7) 162 Dy reaction agrees very well 
with the number of known levels per excitation energy bin below ~1.2 MeV, 
whereas the extracted level density from the 162 Dy( 3 He, 3 He'7) 162 Dy reaction 
overestimates the number of levels in this energy region by a factor of ~3. 
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The level density at ~0.5 MeV of excitation energy is determined by the 
data in the primary 7 matrix which lie approximately on the diagonal E x ~ E^ 
(see Fig. Careful examination of Fig. || shows, that the bumps at E x ~ B 7 
are very well fitted by the factorization given by Eq. ([!]). We therefore con- 
clude, that the differences in level density around ~0.5 MeV of excitation en- 
ergy are not artifacts of the extraction method, but have their origin in dif- 
ferences of the primary 7 spectra. We actually find in the primary 7 ma- 
trix of the 162 Dy( 3 He, 3 He'7) 162 Dy reaction a large number of high energetic 7 
transitions, connecting the direct populated states with the ground state rota- 
tional band. This surplus of counts compared to primary 7 spectra from the 
163 Dy( 3 He,a7) 162 Dy reaction is the reason for overestimating the level density 
at ~0.5 MeV of excitation energy. 

We argue that the level density curve extracted from the neutron pick up 
reaction data is the more realistic one, as supported by Fig. ^. Since the neutron 
pick up reaction cross section is dominated by high I neutron transfer, the 
direct population of the 162 Dy nucleus takes place through one particle one hole 
components of the wave functions. Such configurations are not eigenstates of 
the nucleus, but they are rather distributed over virtually all eigenstates in the 
neighboring excitation energy region. Thus, we can expect fast and complete 
thermalization before 7 emission. The inelastic 3 He scattering on the other hand 
is known to populate mainly collective excitations. These collective excitations 
will thermalize rather slowly, since their structure is much more like eigenstates 
of the nucleus, and their wave functions are less spread over eigenfunctions in the 
close excitation energy region. However, we can expect that their structure is 
similar to the structure of states in the ground state rotational band. Therefore, 
the large 7 transition rates from the direct populated states to the ground state 
rotational band might just reflect the inverse process of inelastic scattering. The 
surplus of 7 counts can therefore be interpreted as preequilibrium decay. An 
extreme example for this are nuclear resonance fluorescence studies (NRF) pl| . 
It is estimated, that in even even nuclei more than 90% of the 7 strength from 
states excited by 7 rays is going to the ground state or to the first excited state. 
Thermalization of the excited states in NRF is also hindered by the fact that 
one populates isovector states, which in the proton neutron interacting boson 
model (IBA-2) are characterized by a different (approximate) F spin quantum 
number than other states in the same excitation energy regions. 

We would like to point out, that although the basic assumption behind the 
primary 7 method is partially violated in the case of the 162 Dy( 3 He, 3 He'7) 162 Dy 
reaction, the level densities in the excitation energy interval 1.5 MeV to 6.5 MeV 
deduced from the two reactions agree extremely well. This indicates, that the 
extracted level density curves are quite robust with respect to the goodness of 
the assumption. Especially the bump at ~2.5 MeV excitation energy indicating 
the breaking of nucleon pairs (?], |l^| and the quenching of pairing correlations 
[|| could be very well reproduced. One should also keep in mind that the 
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two reactions populate states with slightly different spin distributions due to 
the different target spins in the two reactions, which might account for some 
differences in the extracted level densities. 

3.2 The 7 energy dependent factor 

We compare the extracted 7 energy dependent function F of 162 Dy for the two 
reactions. The F function from the 163 Dy( 3 He,a7) 162 Dy reaction was already 
published in p2] using an approximate extraction method, however the data 
were reanalyzed using the exact extraction method of Ref. || and are in the 
present form, as well as data from the 162 Dy( 3 He, 3 He'7) 162 Dy reaction, pub- 
lished for the first time in this work. Figure ^ shows the relative F functions, 
which are obtained by dividing the extracted F function by E™ with n = 4.3 
and scaling them to ~1 at ~4 MeV of 7 energy. Also in this case the two 
functions agree within 10% in the 7 energy interval of 1.5 MeV to 6.5 MeV. 
Above ~6.5 MeV again, the error bars are too large in order to allow for any 
conclusions. Below ~1.3 MeV of 7 energy, the two functions differ dramatically 
from each other. Due to experimental difficulties, like ADC threshold walk and 
bad timing properties of low energetic 7 rays, we had to exclude 7 rays with 
energies below 1 MeV from the data analysis || . It is therefore very difficult to 
judge if the differences in the F function curves below 1.5 MeV of 7 energy are 
also due to experimental problems (i.e. the experimental cut was too optimistic, 
and we should rather have excluded all 7 rays with energies below 1.5 MeV) or 
due to the different nuclear reactions used to excite the 162 Dy nucleus. 

Also here we would like to emphasize, that despite the basic assumption 
behind the primary 7 method is not completely fulfilled in the case of the 
162 Dy( 3 He, 3 He'7) 162 Dy reaction, the two F functions agree very well. Especially 
the bump at ~2.5 MeV of 7 energy, which we interpret as a Pigmy Resonance 
is equally pronounced in both reactions. We are therefore very confident that 
the extracted level density and 7 energy dependent factor for 162 Dy presented 
in this work are not, or very little, reaction dependent. 

4 Conclusions 

This work compares the results from the 162 Dy( 3 He, 3 He'7) 162 Dy reaction to 
those of the 163 Dy( 3 He,a7) 162 Dy reaction. The level density g and the 7 energy 
dependent factor F in 162 Dy are shown to be reliably extracted with our method 
in the energy interval 1.5-6.5 MeV. The findings are independent of the partic- 
ular reaction chosen to excite the 162 Dy nucleus. The two reactions differ from 
each other (i) in the reaction type; i.e. inelastic 3 He scattering versus neutron 
pick up, and thus in the nuclear states populated before thermalization, namely 
collective excitations versus one particle one hole states, (ii) in the target spins; 
+ for 162 Dy versus 5/2~ for 163 Dy, and thus in the spin distribution of direct 
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populated states, and (iii) in the Q-value; MeV for inelastic 3 He scattering 
versus 14.3 MeV for the neutron pick up reaction. Nevertheless, the only differ- 
ences in the extracted quantities are those in the level densities below ~1.5 MeV 
of excitation energy. These might be explained by preequilibrium 7 decay in 
the 162 Dy( 3 He, 3 He'7) 162 Dy reaction, whereas the 163 Dy( 3 He,a7) 162 Dy reaction 
is supposed to show only equilibrium 7 decay, and thus reveals reliable level 
densities below 1.5 MeV of excitation energy, which is supported by comparison 
to known data. However, although preequilibrium 7 decay violates the basic 
assumption of the primary 7 method, the effect on the extracted level density g 
and the 7 energy dependent factor F between 1.5 MeV and 6.5 MeV of energy is 
shown to be less than 10%. In conclusion, the present results have given further 
confidence in the new extraction techniques, and open for several interesting 
applications in the future. 

The preequilibrium decay does not seem to violate the Axel Brink hypothe- 
sis, since the respective parts of the primary 7 spectrum could be fitted within 
this assumption. However, the extracted quantities g and F will then only rep- 
resent a weighted sum of the respective quantities obtained from preequilibrium 
and equilibrium 7 decay, where in the case of the 162 Dy( 3 He, 3 He'7) 162 Dy re- 
action, the preequilibrium process dominates the level density below 1.5 MeV 
of excitation energy. We conclude therefore that neutron pick up reactions are 
more suitable than inelastic 3 He scattering for our method, since the states 
populated by the former reaction presumably thcrmalizc completely, whereas 
those populated by the latter reaction might not completely thermalize before 
7 emission. 
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Figure 1: Raw data of the 162 Dy( 3 He, 3 He'7) 162 Dy reaction. Some 7 spectra for 
different excitation energies are projected out on the right hand side. 




E 7 [MeV] E 7 [MeV] 

Figure 2: Unfolded data of the 162 Dy( 3 He, 3 He'7) 162 Dy reaction. Also here, 
some 7 spectra for different excitation energies are projected out on the right 
hand side. 
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Figure 3: Primary 7 spectra of the 162 Dy( 3 He, 3 He'7) 162 Dy reaction. Here again, 
some 7 spectra for different excitation energies are projected out on the right 
hand side. 
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Figure 4: Normalized primary 7 spectra for the 162 Dy( 3 He, 3 He'7) 162 Dy reaction 
including estimated errors (data points) compared to the least \ 2 fit according 
to Eq. (gj) (lines). 
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Figure 5: Determination of parameters A and a of Eq. (g). The extracted level 
density curve from the 162 Dy( 3 He, 3 He'7) 162 Dy reaction data (full data points 
and line in insert) is compared to the number of known levels per excitation 
energy bin around the ground state (histogram) in the region between the ar- 
rows, and to the level density at the neutron binding energy B n , calculated 
from neutron resonance spacing data (square in insert). In comparison, the ex- 
tracted level density curve from the 163 Dy( 3 He,a7) 162 Dy reaction data (empty 
data points and slashed line in insert) is shown. 
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Figure 6: Comparison of the extracted relative level density of 162 Dy de- 
duced from the 162 Dy( 3 He, 3 He'7) 162 Dy reaction (this work) and from the 
163 Dy( 3 He,cry) 162 Dy reaction (previous works). The error bars of the former 
level density curve are about half of the errors of the latter due to ^5 times 
better statistics in the data of the 162 Dy( 3 He, 3 He'7) 162 Dy reaction. The differ- 
ences between the two curves below 1.5 MeV of excitation energy are accounted 
for in the text. 
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Figure 7: Comparison of the extracted relative 7 energy dependent function 
of 162 Dy deduced from the 162 Dy( 3 He, 3 He'7) 162 Dy (this work) and from the 
163 Dy( 3 He,a7) 162 Dy reaction (previous work, reanalyzed in this work). Also 
here, the error bars of the relative 7 energy dependent function extracted from 
the data of the former reaction are about half of the other ones, due to better 
statistics. 
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